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I. INTRODUCTION 



Gravitational waves emitted by coalescing compact binaries in their adiabatic inspiral phase have been accurately 
modeled [Q and are preferred candidates for the direct detection of gravitational radiation of cosmic origin, by the 
first generation of broadband laser interferometric antennas, including TAMA300 0, GEO600 ||, the LIGOs Q and 
VIRGO §. 

The best strategy for detecting signals of known form in colored gaussian stationary noise consists in correlating 
the detector output (data) with a discrete family (lattice) of expected waveforms (templates), and using the largest 
correlation as a detection statistic (maximum likelihood |6|). The above correlation is the (weighted) scalar product: 



(a, 9) 



c.c. 



(1.1) 



where (/j, / s ) is the antenna spectral window, a(f) — h(f) + n(f) are the (noise-corrupted, spectral) data, h(f) is a 
(possibly null) signal and n(f) a realization of the antenna noise, g(f) is a template, LT(/) is the (one-sided) antenna 
noise power spectral density and C.C. denotes the complex conjugate. The detection statistic has to be compared to 
a threshold, which depends on the allowed false alarm probability. At any fixed false-alarm probability, the largest 
probability of detection is obtained iff a template is exactly matched to the signal. The matching between h and g is 
measured by their overlap, 
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where the norm ||u|| = {u, m) 1//2 . The template lattice should be designed in such a way that for any admissible signal, 
one can always find at least one template in the lattice such that H the overlap is never less than a prescribed T. 

In the straightforward restricted S post-newtonian stationary-phase approximation || the spectral form of a null- 
eccentricity binary inspiral signal is: 



h(f) = Af- 7 ^ cxp [j [27r/T c - C + *M (/, |) - tt/4] } , 



(1.3) 



where A is a (real, unknown) constant, depending on the source position, T c is the fiducial coalescency time |T( 
is the phase at t = T c , and pl| : 
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where v is the PN order, and £ represents the remaining (intrinsic) source parameters. The functions Q and l are 
presently known up to 2.5PN order []l2| and are collected in TableT |D| in terms of the intrinsic parameters m (total 
mass of the binary), rj = fi/m being the reduced mass), f3 and a (spin-orbit and spin-spin terms, respectively 0). 

Maximizing explicitly the overlap over the (irrelevant, extrinsic) parameters A</> c and AT C , one can write the 
template-lattice design prescription as: 
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where the partially maximised overlap M(f,g) is called the match, S and £ are the signal-space and template-lattice, 
respectively, and: 
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Whenever the match is only a function of the signal/template parameter differences, the template lattice is uniform, 
i.e., the lattice spacing is constant throughout the template parameter space, and depends uniquely on V Jl5| . This 
is obviously the case iff the function A^^ depends, in turn, on the signal/template parameter differences only. 

Uniformity is a key property for efficient implementation of the lattice. In a uniform lattice, templates corresponding 
to adjacent nodes differ only by a fixed phase factor, which makes the template family construction computationally 
inexpensive. Furthermore, in view of the quasi band-limited nature of the match as a function of the signal/template 
parameter differences, cardinal interpolation among the correlators can be used as shown in 16 , |fl7| to substantially 
reduce the number of templates required for a prescribed T. 

The ID newtonian (y = 0) and 2D first-order post-newtonian (y = 1) waveforms can be easily parametrized so as 
to make the phase-difference ( |l.6| ) a function of the signal/template parameter differences only 1^]. On the other 
hand, according to present understanding, in order to keep T > 0.97 when working with true data, 2PN (or higher) 
templates will be needed, although spin- free ones should be adequate pl|. 

In the following we shall refer to 2PN templates. These are the best available ones for the purpose of detec- 
tion/estimation. Indeed, as shown in pfj, 2.5PN templates yield generally poorer overlaps (and larger biases) with 
(exact) numerically generated waveforms as compared to 2PN, as an effect of the peculiar (oscillating) behaviour of 
PN-approximant sequences. 

Unfortunately, at PN orders v > 1, as show n in Sect. II, it is strictly impossible to parametrize the spin-free 



waveforms so as to make the phase-difference (L£) a function of the signal/template parameter differences only. 
As a result, the choice of template placement and spacing becomes, in principle, rather unwieldy, and the cardinal 
interpolation techniques discussed in pq| , [PT]] cannot be applied in any straightforward fashion. 

Recently, Tanaka and Tagoshi pl| ] suggested a clever way to circumvent this difficulty, using a simple and elegant 
geometric argument. The scope of this paper is to illustrate the effectiveness of the Tanaka- Tagoshi construction for all 
first generation interferometers (TAMA300, GEO600, LIGO-I and VIRGO) under a given minimal match constraint 
(see Section III). In the following we use geometrized units (G = c = 1) throughout. 



II. THE INTRINSIC CURVATURE OF THE SPIN-FREE PARAMETER SPACE MANIFOLD 

The Tanaka- Tagoshi construction is best understood in the geometric language first introduced in p2| . The match 
between h and g can be accordingly written [^3| : 

M(h,g) = l-G rs A8 r A0 s + ..., (2.1) 



where A8 l = % (^h) — @ l (£,g), and the (2v + l)-dimensional metric G rs is defined in Appendix A. 

Let m m i n < mi < < m max the spin-free waveform parameter space, mi , mi being the companion masses [^4j . 
In the (2iH-l)-dimcnsional PN space 6 l the parameter spa ce is mapped into a three- vertex 2D manifold V . Exploiting 
the dependence of ff' on mi and m-2 one can rewrite (2.1) as: 



M(h, g) = l- ds 2 , ds 2 = g pq Am q Am p , 



(2.2) 
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where 

dO r 09 s 

9pq = Grs dm~ p dm~ q ' (2 - 3) 

is a 2D Finsler metric pf| on V . 

In order to find a chirp waveform parametrization which makes the match a function of the source-template 
'parameter differences only, one should be able to find a coordinate transformation (7711,7712) — > (21,0:2) such that: 

Af (ft, g) = 1 - S pq Ax p Ax q , (2.4) 

where S pq is the 2D euclidean metric. As anticipated, the required coordinate transformation does not exist, in general. 
This is due to the fact that the post-lPN manifold V is not globally flat, as can be seen synthetically from its Gaussian 
curvature K j2(|. The gaussian curvature of V as a function of mi, mi in 0.2M© < mi, 2 < IQM®. is shown in Fig. 1, 
for the special case of LIGOT at 2PN order. The absolute curvature is maximum for mi =7712. This property turns 
out to be common to all first-generation antennas, whose relevant parameters have been collected in Table-II. In Fig. 
2 the curvature of V is displayed as a function of mi =m 2 (worst case) for TAMA300, GEO600, LIGO-I and VIRGO 
at 2PN order. 



III. THE TANAKA-TAGOSHI COORDINATES 



Even though V is not flat, following Tanaka and Tagoshi it is still possible to place the templates on a flat manifold 
chosen close to V in a suitable sense. The Tanaka- Tagoshi construction can be simply phrased as follows: i) introduce 
a linear coordinate transformation 6 — > 2 such that all three vertices of the manifold V are brought onto the (x\, X2) 
plane of the (orthogonal) coordinate system Xi (see Appendix-B for details) , and then ii) take the resulting three- vertex 
flat 2D simplex T in the (21,22) plane as the signal/template parameter space. 

The flat simplexes T corresponding to O.2M < m x < m 2 < 1OM are shown in Fig. 3 for TAMA300, GEO600, 
LIGO-I and VIRGO at 2PN order. Their measures (areas) are collected in Table-Ill. 

Enforcing a minimal-match constraint, in the (21,22) plane yields the uniform square-mesh lattice sidelength 

A= [2(1-T)] 1/2 . (3.1) 

The simplex areas divided by A 2 (area spanned by each template at a given T) provide close estimates of the 
corresponding total number of templates required, also listed in Table-Ill for T = 0.97, for TAMA300, GEO600, 
LIGO-I and VIRGO at 2PN order. 



A. The Minimal Match Error 



Pictorial representations of the departure of V from T are shown Fig. 4, where the simplex T corresponding to 
0.2M Q < 771 1 < 7712 < IOM0 is displayed together with its (euclidean) distance 5 — (xf + x\ + 2 2 ,) 1 / 2 from V as a 
function of 21, x 2 , for TAMA300, GEO600, LIGO-I and VIRGO at 2PN order. The obvious question is to what extent 
does the shown difference spoil the minimal-match condition enforced in T, due to the fact that the true waveform 
space is V. 

Let h a (y— PN) chirp signal, and let H = (x[ h \x 2 h \ . . . ,2^ +1 ) and H = (x[ h \ x^\ 0, . . . , 0) the corresponding 
points in V and T. By construction, V ft. in the admissible spin- free source parameter range m m i n < mi < m 2 < m max , 
there will be at least one lattice node G = (2^ , x^ , 0, . . . , 0) G T such that: 

M(h,g) = 1 - 8 pq {xp-xf){xf)-xf) - T, (3.2) 

where h and g are the waveforms corresponding to H , G, respectively. On the other hand, from (|]l]) , ( |2.3|) 

M(h,g) = l- 7pq ( x W- x M)( x W- x M) = l- n „{x^-xf){xf-xf) (3.3) 

where p,q = l, 2, and p7[ : 

^ = Grs dx^dx^ (3 - 4) 
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where r, s = 1, 2, . . . , 2v + 1. Let x —>■ y the (unit) rotation which diagonalizes the matrix 7 pg — <5 pg . Then: 



M(h,g)-M(h,g) = m (y[ h) -yff+mivT ~vT) 



,{h) „,(5)n2 



< 
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where 771,772 are the eigenvalues of the matrix 7 P9 — <5 P9 . The r-independent quantity: 



77 = 



M(h,g)-M(h,g) 



i-r 



= max(|77i|, I772I) 



(3.6) 



is a measure of the match degradation due to using a flattened parameter space, and is displayed in Fig. 5 as a 
function of mi, ni2 in O.2M < 7711,2 < 10M Q , for the special case of LIGO-I at 2PN order. The maximum value 
is attained for mi =7722. This property is common to all first-generation antennas quoted in Tablc-I. In Fig. 6 the 
quantity 77 is displayed as a function of toi=to 2 (worst case) for TAMA300, GEO600, LIGO-I and VIRGO at 2PN 
order. 

It is seen that the minimal-match condition enforced using T as the waveform/template parameter space is not 
spoiled significantly when the signal (and the template) belongs to P. 

The above argument might perhaps be loosey at low values of T, as e.g. typically required for the initial steps of 
hierarchical search procedures |2S[ ], where the quadratic approximation (2.1) might be no longer accurate enough p3[ . 
In order to evaluate the (maximum) minimal-match error under these broader conditions, we generated 10 4 random 
pairs (777.1,7712), uniformly distributed in the allowed source parameter range 0.2Mq < mi < 7772 < 10M@. For each of 
the above, we drew a closed-curve CcT centered around the corresponding point F E T, whose points G correspond 
to (all) templates g for which M(h,g) = T, and computed 



sup 

Gee 



M(h,g)-M(h,g) 



sup 



1 - 



M(h 7 g) 



(3.7) 



The cumulative distributions of e corresponding to T = 0.7, 0.8, 0.9, for the special case of LIGO-I at 2PN order, are 
shown in Fig. 7. Not unexpectedly, the Tanaka-Tagoshi coordinates are seen to work also at relatively low values of 

r. 



B. Cardinal Interpolation Beyond 1PN 

In the obtained uniform lattice, the application of the cardinal interpolation technique introduced in Jl^ j is straight- 
forward. For 1PN templates this latter yields a fourfold reduction in the template density and total number at T = .97, 
as shown in fl7j| . A main motivation of this study has been to check whether a comparable reduction in the number 
of templates could be obtained using Tanaka-Tagoshi parametrized post-lPN templates. Numerical simulations show 
that this is indeed the case: the 2PN template density reduction as a function of the minimal match T is shown in 
Fig. 8. 



C. Computation of Templates 

In the following we shall denote the (linear) direct (8 — > x) and inverse (x — » 6) Tagoshi-Tanaka transformation 
operators as and 6 U , respectively. 

The template corresponding to the node {x\^ 1 %2,k) of the (uniform) lattice in the (xi, X2) plane should be computed, 
in principle, by solving (e.g., numerically) the system: 

E li 9 l (m lk ,m 2 k) = x\.k ^ ^ 



to obtain the corresponding values mifc,7B2fc of m\,m%, whereby all the Q % functions and hence the (exact) template 
phase (1.4) can be computed. The above exact (and computationally expensive) procedure for computing the template 
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phases needs not to be applied for all lattice nodes. Indeed, as noted in plj] , the change in the d % between the templates 
at {x\,x-i) and (xi + Sxi,X2 + 6x2) is well approximated by: 

A9 l » Q a Sxi + e l2 Sx 2 (3.9) 

as a consequence of the fact that £3,0:4,... ,X2v+i are almost constant throughout relatively large portions of the 
(xi,X2) plane. As a result, the phase difference between neighbouring templates can be taken as uniform throughout 
extended regions of the (xi, X2) plane, which makes the template set construction relatively cheap. One can accordingly 
compute the exact template phases using Eq. (p^q) for a set of sparse lattice nodes. The phase of templates belonging 



to neighbourhoods of these sparse nodes can be obtained using Eq. (3J3). The size of the above neighbourhoods (i.e. 
equivalently, the distance between the sparse nodes) should be chosen so as not to spoil the prescribed minimal match 
condition |29| . 

A possible natural strategy is to place the exact templates nearby the median-lines of T, starting from their 
common crossing point J30|. We checked by numerical simulation that a few (e.g., ~ 50 for LIGO-I at T = 0.97) 



exactly computed sparse templates are sufficient to compute all the remaining ones via (3J3), without spoiling the 
minimal match condition throughout the whole waveform parameter range. 



IV. CONCLUSIONS 



It has been shown that the Tanaka-Tagoshi coordinates are effective to set up a uniform lattice of 2PN (and 
possibly higher) spin-free binary-inspiral templates under a given minimal match constraint, for all 1st generation 
intcrferometric antennas. This allows to extend readily the minimum redundant cardinal interpolation techniques of 
the correlator bank introduced by the Authors to the highest available order PN templates, yielding a computational 
gain comparable to the 1PN case. 

Several possible variations of the Tanaka-Tagoshi method can be envisaged. For instance, a 2D simplex is used here 
to approximate the parameter manifold. Using a 2D complex could possibly improve the accuracy. Further, the 2D 
simplex is obtained by point-collocation, i.e. by enforcing contact between the vertices of the simplex T and those of 
the manifold V. Other approximation philosophies could be used, including e.g., least mean squares, etc. 
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r pq 5x p 5x q < 
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Appendix A - The PN Waveform Parameter Manifold 



In this Appendix we rephrase the well known construction of a metric in the binary-inspiral waveform parameter- 
space, first introduced by Owen [Q and adopted by Tanaka and Tagoshi [^l|, using a slightly different notation. 
Define: 



u(h,g) 



V- 7 / 3 exp{j [27r/AT c -A&+A*M (/;&,&)]} 



df 

n(/) 



c-7/3 



df 

n(/) 



(Al) 
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and let: 



rf. t-7/3 

M/))= '\t. 1-7/3 • ^ 

It is readily proved that: 

u(fc, ff ) = 1 +j{tn)A9 n - i(C m C„)Ar i Ar + . . . , (A3) 
where m, n = 0, 1, . . . , 1v + 1, and 

Co = 2tt(/// ), 6 = f T c . (M) 

From equations (Al) — (A4) one gets: 

max<M> = = 1 - G' mn A9 m A9 n , (A5) 



where m, n = 0, 1, . . . , 2v + 1 and 



2 

Equations (A5) and (A6) allow to express the match as follows: 



G' mn = ~ ((CmCn) - (Cm><Cn» ■ (A6) 



M(ft,p) = maxmax(/i, 5 ) = max |u(/i, 5)] = 1 - G pq A9 p A6 q , (A7) 



where now p,q = 1, 2, . . . , 2u + 1 and pTl 



•"P9 W P9 pi/ 
U 00 



G pg — G' p q . (A8) 



Appendix B - Constructing the Tanaka-Tagoshi Coordinates 

Let m m i n < mi < m,2 < m max the allowed mass range, and let pi, i — 1,2,3 the vertices of the parameter 
space manifold V corresponding to (mi — m2 — m m i n ), (mi=m m i ni m2 = m max ) and (mi—m2 — m max ), respectively, 
in the (2v + l)-dimensional PN parameter space with metric G pq . The Tanaka-Tagoshi coordinates are obtained as 
follows. One first performs a (linear) coordinate transformation 

6i = A} /2 Ve (Bl) 

to make the new coordinates orthogonal. The required transformation is obtained from the Jordan decomposition of 
G, 

G = P T AP, (B2) 

i.e., in component notation: 

G M = P( m) A mn pW, (B3) 

where A mn = X^8 mn , A 1 -™' being the eigenvalues and Pq 11 ^ the eigenvectors of G. 

Next, one seeks a rotation x = Q6*/ which bears (modulo a trivial translation) all vertices p i = 1,2,3 of the 
manifold V onto the coordinate plane (a; 1 , a; 2 ). The problem is thus reduced to that of finding an orthogonal matrix 
Q such that: 

Q(P'3 - P'x) = "nil (B4) 
Q{p f 2 - P'l) = a 2 i£i + a 2 2X 2 , 
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where an, ot2i, a.22 are suitable real numbers, and Xi denote the (new, unit) basis vectors. The system (-B4) does 
not have an unique solution for Q. However a possible solution is readily obtained as follows. Let Z the matrix 
constructed out of the (column) vectors: 

P'3 -P'l, P'2 -P'l, 0/ 3 , ...,§/ 2v +l (B5) 

It is readily shown that Z is non-singular. The straightforward (unique) QR decomposition of Z, viz. [ j32| : 

Z = Q T R, (B6) 

where R is a lower triangular matrix, is obviously a solution of: 

QZ = R. (B7) 

It is seen that the solution of eq. (B7) is also a solution of (BA), and hence the sought Q in {BA) is the same as Q 
in (B6). 
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Table I - Relevant to equation (|l.4j) 
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LIGO-I 


f (^) 4 + 2 




2 


4.4- 10~ 46 


175 


40 


1300 


VIRGO 


n 

4 


290 




+ 1 + 


Jo, 


) 2 




1.1 ■ 10~ 45 


475 


16 


2750 



Table - II - Spectral windows and noise power spectral densities 
of first generation interferometric GW antennas |lq]. 



Antenna 


Simplex Area [sec 2 ] 


No. of Templates at T = 0.97 


TAMA300 


5987 


9.98 x 10 4 


GEO600 


38931 


6.49 x 10 5 


LIGO-I 


16842 


2.81 x 10 5 


VIRGO 


546689 


9.11 x 10 6 



Table III - 2PN flat simplex areas and number of templates at T = .97, 
for O.2M < mi < m 2 < 1OM . 
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CAPTIONS TO THE FIGURES 



Fig. 1 - The gaussian curvature K of V vs. mi, ni2 for LIGO-I. 



Fig. 2 - The gaussian curvature K of V vs. mi = m 2 (worst case) for TAMA300, GEO600, LIGO-I and VIRGO. 



Fig. 3 - The flat simplexes T corresponding to O.2Af < mi < m 2 < 10M Q for TAMA300, GEO600, LIGO-I and 
VIRGO. 



Fig. 4 - The euclidean distance 8 = {x\ + x\ + x^) 1 ^ 2 between the manifolds V and T as a function of x\, X2 for 
TAMA300, GEO600, LIGO-I and VIRGO. 



Fig. 5 - The quantity rj, eq. (3.6), vs. mi, m 2 for LIGO-I. 



Fig. 6 - The quantity r), eq. (EM, vs. mi=m 2 (worst case) for TAMA300, GEO600, LIGO-I and VIRGO. 



Fig. 7 - Cumulative distribution of e, eq. (|j), for LIGO-I at T = 0.7, 0.8, 0.9. 10 4 trial sources in 0.2M Q < mi < 
m 2 < 10M Q . 



Fig. 8 - Template density reduction vs. T after cardinal interpolation (2PN order templates, LIGO-I). 



LIGO-I 




Fig. 1 - The gaussian curvature K of P vs. mi, m 2 for LIGO-I. 



40. 



30. 



20. 



10. 



0. 



-10. 



VIRGO (dots) 
LIGO-I (cont.) 
GEO600 (dots -dash) 
TAMA300 (dash) 




0. 



2. 4. 6. 

mi=m2 |M ] 



8. 



10. 



Fig. 2 - Gaussian curvature K of P vs. mi=m 2 (worst case) 
for TAMA300, GEO600, LIGO-I and VIRGO. 




Fig. 3 - The flat simplexes corresponding to 0.2 M Q < ml < m2 < 10 M 
For TAMA300, GEO600, LIGO-I and VIRGO 




Fig. 4 - The euclidean distance 5 - (x 3 2 + x 4 2 + x 5 2 ) 1/2 between the manifolds P 
and Tzs a function of (x b x 2 ) for TAMA300, GEO600, LIGO-I and VIRGO. 



LIGO-I 




Fig. 5 - The quantity r\ in eq. (III. 6) vs. nil , m 2 for LIGO-I. 



5. 



TAMA300 (dash) 
GEO600 (dots-dash) 
LIGO-I (cont.) 
VIRGO (dots) 




0. 2. 4. 6. 8. 10. 

mi=m 2 [M Q ] 



6 - The quantity r\ in eq. (III. 6) vs. mi = m 2 (worst case), 
for TAMA300, GEO600, LIGO-I and VIRGO. 




Fig. 7 - Cumulative distribution of e , eq. (III.7), for LIGO-I at Y = 0.7, 0.8, 0.9. 
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10 trial sources in 0.2 M < ml < m2 < 10 M 




Fig. 8 - Template density reduction vs. T after cardinal interpolation. 
(2PN order templates, LIGO-I). 



